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Abstract. 

The dynamics of transient electric fields generated by the interaction of high 
intensity laser pulses with undcrdense plasmas has been studied experimentally with 
the proton projection imaging technique. The formation of a charged channel, 
the propagation of its front edge and the late electric field evolution have been 
characterised with high temporal and spatial resolution. Particle-in-cell simulations 
and an electrostatic, ponderomotive model reproduce the experimental features and 
trace them back to the ponderomotive expulsion of electrons and the subsequent ion 
acceleration. 
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1. Introduction 

The study of the propagation of intense laser pulses in underdense plasmas is relevant to 
several highly advanced applications, including electron [1] and ion acceleration [2], [3], 
development of X- and 7-ray sources [4], and fusion neutron production [5]. It is also of 
fundamental interest, due to the variety of relativistic and nonlinear phenomena which 
arise in the laser-plasma interaction [6]. Among these, self-focusing and self-channeling 
of the laser pulse arise in this regime from the intensity dependence of the relativistic 
index of refraction [8] . 

Strong space charge electric fields are generated during the early stage of the 
propagation of a superintense laser pulse through an underdense plasma as the 
ponderomotive force acts on electrons, pushing them away from the axis. Thus, for 
a transient stage the pulse may propagate self-guided in a charged channel [9], while 
the space-charge field in turn drags and accelerates the ions to MeV energies j3]. So far, 
experiments have provided evidence of channel formation and explosion using optical 
diagnostics QUI HH H21 EH E] . or by detecting radially accelerated ions [2j El [H] , 
while a direct detection of the space-charge fields has not been obtained yet. The 
development of the the proton projection imaging (PPI) technique [15] has provided a 
very powerful tool to explore the fast dynamics of plasma phenomena via the detection 
of the associated transient electric field structures. The technique is based on the use 
of laser-accelerated multi-MeV protons ([161 E] an d references therein) as a charged 
probe beam of transient electromagnetic fields in plasmas, a possibility allowed by the 
low emittance and high laminarity of the proton source [HI [19], as well as by its ultra- 
short duration and the straightforward synchronization with an interaction laser pulse. 
The experimental PPI implementation takes advantage from the broad energy spectrum 
of protons, since in a time-of-flight arrangement protons of different energy will probe the 
plasma at different times, and thus an energy-resolved monitoring of the proton probe 
profile allows to obtain single-shot, multi-frame temporal scans of the interaction |15j . 
PPI and the related "proton deflectometry" technique permit to gather spatial and 
temporal maps of the electric fields in the plasma, and therefore have proven to be an 
unique tool to explore the picosecond dynamics of laser-plasma phenomena [201 EH [22] 
via the associated space-charge fields. 

In this article, we report on an experiment using the PPI technique to study of the 
formation and subsequent evolution of a charge-displacement channel in an underdense 
plasma. These investigations have led to the first direct experimental detection of the 
transient electric fields in the channel, providing an insight of the fundamental physical 
processes involved. The comparison of the experimental data with two-dimensional (2D) 
electromagnetic (EM) particle-in-cell (PIC) simulations and a simple one- dimensional 
(ID) electrostatic (ES) PIC model allows to characterize in detail the electric field 
dynamics at different stages of its evolution. 
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2. Experimental setup 




Figure 1. (color) (a) Schematic of the experimental setup, (b) Top: Interferogram 
of the plasma at 25 ps before the arrival of CPAi at its focal plane x — 0. The CPAi 
peak intensity was 1.5 x 10 19 W cm~ 2 . Bottom: the corresponding electron density 
profile along the y — axis. 

The experiment was carried out at the Rutherford Appleton Laboratory, employing 
the VULCAN Nd-Glass laser system [23], providing two Chirped Pulse Amplified (CPA) 
pulses, with 1.054 /im wavelength, synchronized with picosecond precision. Each of the 
beams delivered approximately 30 J on target in 1.3 ps (FWHM) duration. By using 
f/6 off-axis parabolas, the beams were focused to spots of 10 /im (FWHM) achieving 
peak intensities up to 3 x 10 19 W cm -2 . The short pulses were preceded by an Amplified 
Spontaneous Emission (ASE) pedestal of 300 ps duration and contrast ratio of ~ 10 6 
|24j . One of the beams (CPAi) was directed to propagate through He gas from a 
supersonic nozzle, having a 2 mm aperture, driven at 50 bar pressure. The interaction 
was transversely probed by the proton beam produced from the interaction of the second 
CPA beam (CPA2) with a flat foil (a 10 /im thick Au foil was typically used), under the 
point projection imaging scheme [15]. The schematic of the experimental setup is shown 
in FigJTTa). Due to the Bragg peak energy deposition properties of the protons, the 
use of multilayered stacks of Radiochromic film (RCF) detector permits energy-resolved 
monitoring of the proton probe profile, as each layer will primarily detect protons within 
a given energy range. This allows to obtain single-shot, multi-frame temporal scans of 
the interaction in a time-of-flight arrangement [15]. The spatial and temporal resolution 
of each frame were of the order of a few ps, and of a few microns, respectively, while 
the magnification was 11. 

The interaction region was also diagnosed by Nomarsky interferometry, employing 
a frequency doubled CPA pulse of low energy. The reconstructed electron density profile 
along the propagation axis before the high-intensity interaction [see Fig. d^b)], broadly 
consistent with the neutral density profile of the gas jet [25] suggests complete ionisation 
of the gas by the ASE prepulse. 
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3. Experimental Results 

Fig. 12(a), (b) and (c) show three sequential proton-projection images of the interaction 
region. The laser pulse propagates from left to right. A 'white' channel with 'dark' 
boundaries is visible at early times [Figj2](a) and (b)]. The leading, 'bullet' shaped edge 
of the channel, indicated by the label I in Fig [5] (a) and (b), is seen moving along the 
laser axis. In the trail of the channel the proton flux distribution changes qualitatively 
[Fig. [2(c) and (d)], showing a dark line along the axis (indicated by the label III), which 
is observed up to tens of ps after the transit of the peak of the pulse. The 'white' 
channel reveals the presence of a positively charged region around the laser axis, where 
the electric field points outwards. This can be interpreted as the result of the expulsion 
of electrons from the central region. The central dark line observed at later times in 
the channel suggests that at this stage the radial electric field must change its sign, i.e. 
point inwards, at some radial position. As discussed below, this field inversion is related 
to the effects of ion motion. 

Due to multi-frame capability of the PPI technique with ps temporal resolution, 
it has been possible to estimate the propagation velocity v of the channel front. 
Critically, one has to take into account that, due to the divergence of the probe beam, 
the probing time varies along the pulse propagation axis (see Fig. [3]) as t(x, E) ~ 
t (E) + r (E)( v /l + x 2 /L 2 - 1) where L ~ 0.3 cm is the distance between the plane 
and the proton source and r = L / y^2E/m p ~ 220 ps/ ^E/MeV is the proton time 
of flight from the source to the center of the object plane (see Figj3j). The reference 
time to(E) is relative to the instant at which the laser pulse peak crosses the focal 
plane x — 0. We divide the displacement of the tip of the channel leading front from 
frame (a) to (b), AX = X b - X a ~ (-300 + 500) /mi = 200 /xm by the difference At 
in the corresponding probing times At = t(X&, E b ) — r(X a , E a ) ~ 0.66 ps to obtain 
v = AX/ At ~3x 10 s m/s. Based on the nominal reference time we estimate the peak 
of the pulse to be approximately 0.6 ps behind the tip of the channel front. 

4. Data analysis by numerical simulations 

In order to unfold the physical mechanisms associated with the dynamics of the charged 
channel, a ID ES PIC model in cylindrical geometry was employed, in which the laser 
action is modeled solely via the ponderomotive force (PF) of a non-evolving laser pulse. 
A similar approach has been previously used by other authors [T2l [14"] . The code solves 
the equation of motion for plasma particles along the radial direction, taking into account 
the ES field obtained from Poisson's equation and the PF acting on the electrons - 
F r = -m e c 2 d r [l +a 2 (r,t)/2] 1 / 2 [26]. Here a(r,t) = a e~ r2 / 2r ° f (t) and f(t) defines the 
temporal envelope of the laser pulse. For the latter, a 'sin 2 ' profile was used (the use of 
a Gaussian profile did not yield significant differences). 

Fig. [4] shows the electric field E r (r,t) and the ion density rii(r,t) obtained from a 
simulation with ro = 7.5A, do = 2.7 and initial density 0.01n c . The initial depletion of 




Figure 2. Proton projection images of the interaction region at different times 
obtained in two different laser shots. The x and y coordinates refer to the object 
(interaction) plane, which intersects the probe axis at (x,y) = (0,0). The images (a) 
and (b) have been obtained from the same shot at an intensity / = 4.0 x 10 18 W cm~ 2 , 
while (c) and (d) correspond to a shot with I = 1.5 x 10 19 W cm~ 2 . The signal in 
the frames is mainly due to protons [of energies E = 13 MeV, 12.5 MeV, 13 McV and 
12 MeV in (a)-(d) respectively] reaching the Bragg peak within their active layers. 
The time labels give the probing time of the protons propagating along the probe axis 
[to (E)] , relative to the time of arrival of the peak of the interaction pulse at the plane 
x = /im. White (dark) regions correspond to lower (higher) proton flux than the 
background. The labels indicate the most prominent features: I) the bullet-shaped 
leading edge and II) the central region of the 'white', positively charged channel; III) 
the 'black' line along the axis, indicating a region of field inversion inside the channel. 



electrons and the later formation of an ambipolar electric field front are clearly evident. 
In order to achieve a direct comparison with the experimental data, a 3D particle tracing 
simulation, employing the PTRACE code [22], was carried out to obtain the proton 
images for an electric field given by E(x,r, t) = rE r (r,t — x/c). The experimental 
proton source characteristics (e.g. spectrum and divergence), the detector configuration 
and dose response were taken into account. A simulated proton image, reproducing 
well the main features observed in the experiment, is shown in Fig. H](c). The tip of 
the channel front is located 0.75 ps ahead of the pulse peak, in fair agreement with the 
previous estimate based on the data. 
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Figure 3. Schematic showing the evaluation of the channel front speed. Protons of 
energy E from the point-like source cross the channel axis at a distance x from the 
centre of the object plane at the time t(x,E) = t (E) + [L(x) — L ]/v p , where L is 
the distance between the source and the object plane, v p = \J2E /m p is the proton 
velocity, L(x) = \Jx 2 + Lq, and to(E) is the time at which the protons of energy E 
directed perpendicularly to the object plane cross the latter, relative to the instant at 
which the laser pulse peak reaches the focal plane (x = 0). 

An essential description theoretical description of the dynamics observed in the ID 
simulations can be given as follows (a detailed description is reported in Ref. [27]). In the 
first stage, F r pushes part of the electrons outwards, quickly creating a positively charged 
channel along the axis and a radial ES field which holds the electrons back, balancing 
almost exactly the PF, i.e. eE r ~ F r ; thus, the electrons are in a quasi-equilibrium state, 
and no significant electron heating occurs. Meanwhile, the force ZeE r ~ ZF r accelerates 
the ions producing a depression in m around the axis [see FigJH (b)]; at the end of the 
pulse (t ~ 2 ps) we find E r ~ [see FigJH (a)], indicating that the ions have reached the 
electrons and restored the local charge neutrality. However, the ions retain the velocity 
acquired during the acceleration stage. For r > r max , where r max is the position of the 
PF maximum, the force on the ions, and thus the ion final velocity, decrease with r; 
as a consequence, the ions starting at a position rj(0) > r max are ballistically focused 
towards a narrow region at the edge of the intensity profile, and reach approximately 
the same position (r ~ 15 jum) at the same time (t ~ 3 ps, i.e. 2 ps after the peak of the 
laser pulse); here they pile up producing a very sharp peak of 7ij. Correspondingly, the 
ion phase space shows that the fastest ions (with energies ~ 400 keV, of the order of the 
time-averaged ponderomotive potential) overturn the slowest ones and hydrodynamical 
breaking of the ion fluid occurs. Using a simple model [27], the "breaking" time and 
position can be estimated to be n ~ t v + (n/2\ /r 2)e 3 ^^y(A/Z)(m p /m e )(r /a c) (where 
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Figure 4. Top: simulated proton image, obtained from particle tracing simulations 
of 12 MeV protons [as in the data in Fig. [2] (d)] in the electric field pattern E(r, z, t) 
given by the ID particle simulations based on the ponderomotive, electrostatic model. 
Bottom : profiles of electric field E r and ion density rii, and phase space distributions 
of ions fi(r,p r ) and electrons f e {r,p r ) from the ID simulation, at various times (t = 
refers to the laser pulse peak). The simulation parameters were ao = 2.7, ro = 7.5 fj,m, 
and initial electron density n e = lx 10 19 cm~ 3 . 



t p = 1 ps is the time at which the pulse has maximum amplitude) and r& ~ (3/2) 3 / 2 ro, 
yielding Ty, — t p ~ 1.3 ps and r& ~ 14 /xm for the simulation in FigJU in good agreement 
with the numerical results. As inferred from the ion phase space plot at t = 4.7 ps, 
a few ions acquire negative velocity after "breaking"; they return toward the axis and 
lead to the formation of a local density maximum at r = after t ~ 15 ps. 

The electron phase space shows that at breaking the electrons are strongly 
heated around the ion density peak, generating an "hot" electron population with a 
"temperature" ~ 13 keV and a density ~ 4 x 10 19 cm -3 , corresponding to a local 
Debye length is Ad = (7)i/47rn/ 1 e 2 ) 1 / 2 ~ 0.13 fxm. A modeling of the sheath field thus 
generated around the density spike [27] (whose thickness d ~ 0.1 fim is less than both 
\d and the sheath width L) yields a peak field E s ~ 2nenhd ~ 6 x 10 10 V m _1 and 
a sheath width L ~ 4X 2 D /d ~ 0.7 /iin, consistently with the simulation results. The 
ambipolar field at the "breaking" location can be thus be interpreted as the sheath field 
resulting from the local electron heating. 

The good agreement of the simulated image with the experimental ones indicates 
that the ID ponderomotive, electrostatic model contains the essential physics of self- 
channeling and related electric field dynamics, despite the exclusion of nonlinear pulse 
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Figure 5. Top frame: ion density (nj) and electric field components (E z and E y ) 
at t = 2.0 ps. Bottom frame: lineout of E y (blue) and rii (red) along the y-axis 
at two different x-positions. showing the transition in the radial field pattern. The 
density is normalized to n c = 10 21 cm" 3 and the fields to m e uic/e = 3.2 x fO 11 V/m. 
The initial density profile reproduces the experimental one in FigJT] (b). The laser 
pulse propagates from left to right along the x axis and has 1 ps duration. The peak 
amplitude is ao = 2 in dimensionless units. The laser pulse is s-polarized (i.e. the 
polarization is perpendicular to the simulation plane). In this configuration E z is 
representative of the amplitude of the propagating EM pulse and E y is generated by 
the space-charge displacement. The observed increase in E z by a factor of ~ f .2 with 
respect to the peak value in vacuum is due to self-focusing. 



evolution due, e.g., to self- focusing, or of the plasma inhomogeneity. To address these 
issues we performed electromagnetic particle-in-cell simulations of the laser-plasma 
interaction. The simulation were two-dimensional (2D) in planar geometry - a fully 3D 
simulation with spatial and temporal scales close to the experimental ones and adequate 
numerical resolutions is way beyond present days computational power. The 2D results 
can be considered to be qualitative since, for instance, the diffraction length of the laser 
beam or the scaling with distance of the electrostatic field are different in 3D. Despite 
of these limitations, the main features of the channel observed in the experimental data 
are qualitatively reproduced in 2D simulations for a range of parameters close to the 
experiment. In the simulation of Fig. [5j the laser pulse has a Gaussian intensity profile 
both in space and time, with peak dimensionless amplitude ao = 2, radius ro = 4A 
and duration r = 300 A/c where A is the laser wavelength. For A = 1 «m the pulse 
duration and intensity correspond to 1 ps and 5.5 x 10 18 W cm -2 , respectively. The 
charge-to- mass ratio of ions is ZjA = 1/2. The electron density grows linearly along 
the x-axis from zero to the peak value no = 0.1n c (where n c is the critical density and 
n c = 10 21 cm" 3 for A = 1 urn) over a length of 400A, and then remains uniform for 200A. 
A 6500 x 1200 numerical grid, a spatial resolution Ax = Ay = A/10 and 16 particles 
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per cell for both electrons and ions were used. 

Fig. [5] shows the ion density (n$) and the components E y and E z of the electric 
field at the time t = 600A/c ~ 2.0 ps. In this simulation the laser pulse is s-polarized, 
i.e. the polarization is along the z axis, perpendicular to the simulation plane. Thus, 
in Fig. [5] E z is representative of the amplitude of the propagating EM pulse, while E y 
is generated by the space-charge displacement. Simulations performed for the case of 
p-polarization showed no substantial differences in the electrostatic field pattern. The 
simulation clearly shows the formation of an electron-depleted channel, resulting in an 
outwardly directed radial space-charge electric field whose peak value is 6.7 x 10 10 V/m 
(see the lineout at x — 410A in Figj5]). In the region behind the peak of the pulse, two 
narrow ambipolar fronts (one on either side of the propagation axis) are observed. The 
ambipolar fields have peak values of ~ 6 x 10 10 V/m (see the lineout at x = 325A in 
Figj5]). As shown above, such a radial electric field profile produces a pattern in the 
proton images similar to that observed in region (III) of Fig. [2j 

5. Conclusion 

We have reported the first direct experimental study of the electric field dynamics in a 
charge-displacement channel produced by the interaction of a high intensity laser pulse 
with an underdense plasma. The field profiles observed clearly identify different stages 
of the channel evolution: the electron depletion near the axis due to the ponderomotive 
force, and the following ion acceleration causing a field inversion along the radius. The 
features observed are reproduced and interpreted by means of ID electrostatic and 2D 
electromagnetic PIC simulations, followed by a reconstruction of the proton images 
employing a 3D particle tracing code. 

Acknowledgements 

This work has been supported by an EPSRC grant, Royal Society Joint Project and 
Short Visit Grants, British-Council-MURST-CRUI, TR18 and GRK1203 networks, and 
MIUR (Italy) via a PRIN project. Part of the simulations were performed at CINECA 
(Bologna, Italy) sponsored by the INFM super-computing initiative. We acknowledge 
useful discussions with F. Cornolti and F. Pegoraro and the support of J. Fuchs and the 
staffs at the Central Laser Facility, RAL (UK). 

References 

[1] Malka V. et al, Plasma Phys. Control. Fusion , 47, B481 (2005) and the references therein. 

[2] Krushelnick K. et al, Phys. Rev. Lett. , 83, 737 (1999); Wei M. S. et al., ibid. 93, 155003 (2004). 

[3] Willingale L. et al, Phys. Rev. Lett. , 96, 245002 (2006). 

[4] Roussc A. et al, Phys. Rev. Lett. , 93, 135005 (2004). 

[5] Fritzlcr S. et al, Phys. Rev. Lett. , 89, 165004 (2002). 

[6] Bulanov S. V. et al, Rev. Plasma Phys., 22, 227 (2001). 

[7] Sun G. Z. et al, Phys. Fluids, 30, 526 (1987). 



Dynamics of charge- displacement channeling in intense laser-plasma interactions 10 



Mori W. B. et al, Phys. Rev. Lett. , 60, 1298 (1988). 

Borisov A. B. et al, Phys. Rev. Lett. , 68, 2309 (1992); Borisov A. B. et al, Phys. Rev. A 45, 
5830 (1992). 

Monot P. et al, Phys. Rev. Lett. , 74, 2953 (1995). 
Borghesi M. et al, Phys. Rev. Lett. , 78, 879 (1997). 
Krushelnick K. et al, Phys. Rev. Lett. , 78, 4047 (1997). 
Sarkisov G.S. et al, JETP Lett, 66, 828 (1997). 
Sarkisov G. S. et al, Phys. Rev. E, 59, 7042 (1999). 
Borghesi M. et al, Rev. Sci. Instrum. , 74, 1688 (2003). 
Fuchs J. et al, Nature Phys., 2, 48-54 (2006). 
Robson L. et. al, Nature Phys., 3, 58 (2006). 
Borghesi M. et al, Phys. Rev. Lett. , 92, 055003 (2006). 
Cowan T. et al, Phys. Rev. Lett. , 92, 204801 (2006). 
Borghesi M. et al, Phys. Rev. Lett. , 88, 135002 (2002). 
Borghesi M. et al, Phys. Rev. Lett. , 94, 195003 (2005). 
Romagnani L. et al, Phys. Rev. Lett. , 95, 195001 (2005). 
Danson C. et al, J.Mod.Opt., 45, 1653 (1998). 

Gregori G., Central Laser Facility, RAL, UK, private communication. 

Jung R. et al, in: Central Laser Facility Annual Report 2004/2005, RAL Report No. RAL-TR- 
2005-025, p.23. 

Bauer D., Mulser P. and Stceb W. H., Phys. Rev. Lett. , 94, 195003 (2005). 
Macchi A. et al, physics/0701139. 



